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Abstract 

This is the first part of a survey whose ultimate purpose is to 
clarify the significance of the famous coincidence between the Hubble age 
of the universe and a certain combination of microphysical parameters. In 
this part the way is prepared by a discussion of the manner in which familiar 
local phenomena depend qualitatively, and in order of magnitude, 
quantitatively on the fundamental parameters of microphysics. In order 
to keep the account concise while remaining self contained, only the 
barest essentials of the standard nuclear physical and astrophysical calculat- 
ions involved are given. Only six of the fundamental parameters play a 
dominant part, namely the coupling constants of the strong, electromagnetic, 
and gravitational forces, and the mass ratios of the proton, neutron, electron 
and TT-meson. Attention is drawn to the important consequences of three 
coincidental relationships between these parameters. It is shown that 
most of the principle limiting masses of astrophysics arise (in fundamental 
units) simply as the reciprocal of the gravitational fine structure constant, 
with relatively small adjustment factors. The dividing point between red 
dwarf and blue giant stars turns out to be an exception: this division 
occurs within the range of the main sequence stars only as a consequence of 
the rather exotic coincidence that the ninth power of the electromagnetic 
fine structure constant is roughly equal to the square root of the 
gravitational fine structure constant. 
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1. Introduction 

The main aim of Part I of this survey is to assess the extent to 
which famihar large scale (but local) natural phenomena depend for their 
qualitative and qualitative character on the values of fundamental 
microphysical parameters, and particularly to distinguish those features 
which do not depend critically on the values of these parameters from 
those which depend on numerical coincidences. A secondary aim is to 
prepare for Part II where the discussion will be extended to cosmology with 
the ultimate objective of throwing light on the well known coincidence^ 
which is most commonly expressed as a comparison between the ratio of the 
Hubble radius of the universe to the classical radius of the electron and 
the ratio of the electric and gravitational attraction between a proton 
and an electron. These ratios are both extremely large - of the order 10^" - yet 
they are observed to agree within a factor of order ten. Part II consists 
of an attempt to show that this coincidence can be fully explained in 
principle (although many relevant details remain uncalculated in practice) 
in terms of conventional physics and cosmology, so that revolutionary 
departures such as Dirac's hypothesis of varying gravitational constant,^ 
or Eddington's Fundamental Theory^ are not justified. The connection 
between local and cosmological quantities will be achieved via the timescales 
of stellar evolution, and therefore the final task of Part I is to obtain 
formulae for these timescales in terms of fundamental microphysical 
parameters. 
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In a discussion of this sort it is desirable to express all quantities 
in terms of fundamental units so as to avoid being distracted by dimensional 
numbers which have no direct physical significance. For this reason we 
shall throughout use basic units in which Newton's gravitational constant, G 
the Dirac form h of Planck's constant, and the speed of light, c, are 
set equal to unity. For charges we shall use the corresponding unrationalised 
units, and for temperature we shall use the corresponding units with 
Boltzmann's constant k set equal to unity. We shall deal mainly in 
orders of magnitude, and for this purpose the rough equality symbol ~ 
will be used, implying that the left and right hand sides do not differ by 
much more than a factor of order ten. When greater precision is desired, 
the approximate equality symbol = will be used implying that the 
left and right hand sides do not differ by much more than a fraction of 
order one tenth. 

Throughout the discussion a primary role will be played by a single very 
small number namely the gravitational coupling constant of the nucleon. The 
nucleon mass is = 8 x 10~^° in fundamental units (the 
difference between the proton and neutron masses being negligible in this 
context) which gives rise to a gravitational coupling constant m"^ = 3/5 x 10~' 
this being what is customarily known as the gravitational fine structure 
constant. Nearly all of the very large numbers of cosmogony arise 
essentially as simple small powers of this number. One exceptional case 
will however be found, where a very large number (determining the 
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dividing point between stars with convective and radiative envelopes) 
arises from a high power of a relatively moderate number, namely 
the ordinary (electromagnetic) fine structure constant, = 1/137 where 
e is the electron charge. 

It will be shown that most important limiting masses in astrophysics 
are closely related to the Landau mass. Ml, defined as the reciprocal 
of the gravitational fine structure in fundamental units, i. e. 

Ml = m^' (1) 

It is so named because its importance was first realised by Landau^'^ who 
showed in 1932 (prior to more accurate calculations by Chandrasekhar,^'*^ 
Oppenheimer and Volkhoff^ etc.) that it gives the order of magnitude of the 
largest mass which can support itself as a cold spherical body against 
gravitational collapse. It will appear that other limiting masses than 
this differ from Ml only by a moderate factor whose origin is either 
arithmetical or else dependent on moderate valued fundamental parameters 
such as the electromagnetic fine structure constant. 

This simple and important fact does not seem to have been widely 
appreciated except by specialists in high energy astrophysics. The sort of 
misconception which has existed is illustrated by the theory of Jordan (1947)^. 
Jordan noticed that the number iVj defined as the nucleon content of 
the largest commonly observed types of star, is of order 10^^, and he 
was much impressed by the fact that this is of the same order as the three 
halves power of the ratio (of the Hubble radius to the classical electron 
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radius) which appears in the famous cosmological coincidence. He therefore 

concluded that this would have to be explained by the existence of a 

mechanism whereby matter was originally created directly in the form of 

stars, with a cosmologically determined size limit. However Bondi (1960) 

has pointed out that most astrophysicists believe that the upper limit 

is adequately accounted for by conventional theories of stellar structure, 

although without making clear whether such theories give a direct explanation 

of the coincidence or whether it is accidental. Actually Jordan's coincidence 

is indeed a direct consequence of the cosmological coincidence, since the 

ratio of electromagnetic to gravitational attraction between a proton and 

an electron is e^/rriem-^ , where rrie is the electron mass, 

Jordan's coincidence can be recast in the form A'^j ^ (e^/mNme)^/^ . 

Now this result can only be puzzling if one fails to realise that Nj must 

be closely related to the Landau nucleon number corresponding to the Landau 

mass defined by A'l = Mi^/m^ i.e. by 

Nl = m^' (2) 

Using this coincidence we can rewrite Jordan's coincidence as A^j ~ {e^m^ / rricY/'^ Ni^. 

It will be shown that according to conventional theories of stellar 
structure no normal stable star can exist with a nucleon number which differs 
from the Landau number by more than a factor of order 10^ either way, a 
result which, it should be emphasised, does not depend on how the star was 
formed (which is fortunate because the theory of stellar formation is not 
yet in a very satisfactory state). The lower limit lies so close to A^l 
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only as a consequence of the microphysical coincidence that mg ~ e^m^ 

where is the mass of the vr-meson. However the upper hmit 

hes close to A'l automatically, and in this case the factor 10^ is 

of purely arithmetical origin, so that the physically correct formula for 

Jordan's number is A'^j ^ IO^A^l- Thus it appears that 

Jordan was correct in guessing that A'^j is basically the three halves 

power of the large number mj^^ ~ i x 10'^^, although quite 

mistaken in guessing the reason for this dependence. The small factors 

and and ttin / rric turn out to be quite irrelevant in this case - it 

is purely accidental that the combination [e^m-^/mcY^'^ = 16 

is of the same order of magnitude as the rather generous adjustment factor 

102. 

Nearly all the relationships and coincidences which will be mentioned 
in the following sections are well known and understood by experts in the 
topics concerned. They are presented here in a form which is intended to 
be readily accessible and comprehensible to physicists generally, but it 
is also hoped that specialists may be able to glean some new insights from 
the unified treatment which will be given. With this end in view we shall 
go as rapidly as possible through the steps of the relevant order of 
magnitude calculations, picking out the physically dominant quantities at 
each stage. The crudest possible approximations will be used throughout; 
no apology is offered for this since the purpose is to focus attention on 
the physical mechanism most fundamentally involved and to bring to light 
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relationships which are usually obscured by the mass of details which 
result from more sophisticated calculations. (The collection of 
formulas which results can also be used as a guide to the effects which 
can be expected in cosmological theories where the fundamental 
microphysical parameters are functions of time or position, the most commonly 
considered kinds involving variation of Newton's "constant" (which in 
fundamental units means variation of absolute elementary particle masses 
while keeping the mass ratios fixed) and variation of the fine structure 
parameter. Since at each stage detailed theoretical or (in the few 
cases where the theory is unreliable) experimental evidence is available 
as a check, the results can be regarded as sound and reliable except in 
a few cases which will be pointed out as they arise. 

2. Fundamental Microphysical Parameters 

At the present moment there is a rather extensive set of micro- 
physical parameters which are fundamental in the sense that the reasons 
for their values and inter-relationships are beyond the scope of 
currently accepted physical theories. The most important of these, and 
the only ones we shall need to consider here in relation to natural 
phenomena are coupling constants: more specifically the coupling 
constants of the strong, weak, and electromagnetic interactions, and the 
elementary particle masses, which can be thought of as coupling constants 
of the gravitational interaction. 
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In addition to these coupling constants there are many other kinds 
of fundamental microphysical parameters, (for example magnetic moments of 
elementary particles), but because these other kinds do not seem to play 
an important part in macrophysical phenomena we shall not need to take 
account of them. On the other hand it is no longer necessary to treat 
all the elementary particle masses as fundamental, because SU(3) and related 
theories have achieved considerable success in relating mass splittings 
within multiplets^". However there is at present no generally accepted 
method of calculating the absolute values of strong and electromagnetic 
mass splittings (although attempts are being made)^^, and calculations of the 
original masses about which splittings take place are even further 
from being achieved. 

In order to determine the general character of everyday natural 
phenomena (both terrestrial and astronomical) it is only necessary to 
know the parameters connected with the strong electromagnetic, and 
gravitational interactions. This is not to deny that weak interactions 
are, in themselves, of the greatest importance: in fact they play a key 
role in making possible many crucial nuclear reactions. However the 
quantities and rates of energy production will normally be determined by 
parameters arising from other types of interaction. Thus provided we 
know that the neutrino mass and the weak coupling constant are sufficiently 
small (and that the latter is not too small) their exact values do not 
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matter very much, of course these remarks do not apply to the 
cataclysmic processes that may occur in a supernova or in gravitational 
collapse. 

Nearly all the natural terrestrial phenomena (except radioactive decay 
and cosmic ray interactions) depend essentially only on three fundamental 
parameters, the electron charge e, the nucleon mass ttin, and 
the electron mass me, provided that the purely integral values of 
the nucleon numbers and charge numbers of the stable nuclei are taken 
for granted. Additional information, concerning exact atomic weights 
for example would only be necessary for calculating fine corrections and 
would not have any qualitative effect (unless perhaps some biological 
phenomena depends in a subtle manner on very exact energy values - which 
seems rather unlikely). 

However if we wish to calculate the nucleon and charge numbers of 
atoms (not to mention corrections to atomic weights) and if we wish to 
understand common astrophysical phenomena, (or the historical origin of 
the earth), for which thermonuclear reactions are of central importance, 
then we need three more fundamental parameters. These are the pseudoscalar 
coupling constant, gs, of the strong interactions, the pion mass m.„ 
which determines the maximum effective range of the strong interactions, 
and thirdly the electromagnetic mass splitting An, of the nucleon, 
defined as the excess of the neutron mass over the proton mass. 

These six parameters are sufficient to determine the character of all 
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important natural phenomena except in those rare and exotic cases when 
high energy physics is required or when cosmological quantities are 
directly involved. They can be categorised as three coupling constants, 
and three mass ratios, and their empirically determined numerical 
values are approximately:^^ 

4 e = ^ mN = ^x 10-10 (3) 

7 mN ~ 1830 mN ~ 730 ^ ' 

The values of the coupling constants are rather more familiar in their 
squared forms: thus we have the gravitational fine structure constant 

= 2 X 10"'^^ the ordinary (electromagnetic) fine structure 
constant, = 1/137, and what may be called the course structure 
constant, = 15. It is customary to use the absolute mass 
of the nucleon (rather than say the electron or the pion) for defining 
the gravitational coupling constant (and hence the corresponding fine 
structure constant) not because there is any known theoretical reason 
why the nucleon mass should be more fundamental than the others, but 
simply because the masses of most natural objects can be most conveniently 
be accounted for approximately in terms of their nucleon content. (Electrons 
give rise to very small corrections, and pions and other particles exist 
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mostly in virtual states, and give rise to corrections which are seldom 
much larger). 

None of the six numbers differs outrageously from unity except for 
the gravitational coupling constant with value about 10^^^. The main 
theme of this paper will be a demonstration of the way in which most of the 
large numbers of cosmogony arise naturally as powers essentially of this 
number only, the other fundamental parameters giving only relatively fine 
corrections. 

One can complete a list of the more traditional microphysical parameters 
(i.e. the masses and coupling constants known before the discovery of the 
strange particles of the SU(3) multiplets) by adding three more: the 
muon mass, m^, the electromagnetic mass splitting, A^, of the 
pion (defined as the mass excess of the charged over the uncharged pion) 
and the weak coupling constant, (^W; which characterises /3-decay, 
using the nucleon Compton wavelength as a length scale(it is not uncommon 
for the pion Compton radius to be used as a lengthscale giving rise to the 
modified value (m7r/mN)^5'W; but the choice is arbitrary, and 
it is consistent with the scheme adopted here to use the nucleon mass for 
calibration purposes throughout). Their values are given by 

m,, 1 1 c 

— = - — = gw = 10"^ 5 

mN 9 rrtN 200 ^ ^ 

but for the reasons already indicated, these numbers do not have any macro- 
physical importance under normal conditions. 
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Three remarkable relations between the six important parameters 
stand out immediately from the above tabulation: to a good approximation, 
we have (in order of decreasing accuracy), 

= 2me/m^ (6) 
gs = 2m^/m^ (7) 
An = 2me (8) 

These three relations are not only interesting in themselves. We shall 
see in the next and subsequent sections that each one of them is of consider- 
able phenomenological importance. In fact (4) and (5) play such 
a critical role that even the removal of the factor 2 from their right hand 
side would have drastic consequences. 

For the sake of completeness we mention three similar relations 
involving the other three parameters, even though they will not play any 
part in the following sections. They are (again in decreasing order of 
accuracy 

tt^nAn = 2m7rA^ 9) 
~ (10) 
g^^ = 4mN (11) 
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The last of these is the most interesting since it connects the only two 
parameters on the list which differ greatly from unity. It states that 
the square of the weak fine structure constant is of the same magnitude 
as the square root of the gravitational fine structure constant. 

As far as present day microphysical theory is concerned, all these 
relationships are purely fortuitous, although it is naturally to be hoped 
that future developments will be able to explain them. Some people have 
attempted to account for them in terms of cosmological concepts (e.g. 
Hayakawa^^ has offered a heuristic and speculative explanation for the 
relation g^/e^ = m-f^/m^, which arises as the quotient of (7) 
and (8). However most physicists would probably expect that since the 
quantities involved are purely microphysical their explanation should be 
microphysical also. (If any one of them does have a cosmological explan- 
ation it is most likely to be (11) since only this one directly involves 
gravitation). In any case it is by now means unquestionable that they 
have any direct significance at al: they may arise in a roundabout way 
from the interaction of many different effects. Some coincidences of this 
sort may be expected to turn up in any set of arbitrary numbers. Indeed 
one could easily write down more, e.g. by correct choice of the index 

e^^~mN (12) 

(a relation which, rather surprisingly, will be found to have phenomenol- 
ogical importance resulting from the precise value of the index required), 
of course an appearance of precision can always be given by insertion of 
factors such as 7r/2 etc. in the equations, but ad hoc adjustments add little 
to their real significance. Probably the first one to be explained will be 
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(9) since active attempts to calculate electromagnetic mass splittings 
are already being made^^. It might be hoped that if such an attempt were 
successful, an explanation of (5) would soon follow, but this would 
not really be likely since on one hand electrons are not thought to be 
directly involved in the mass splitting effect, while on the other hand 
a calculation of the absolute electron mass is not even on the horizon. 
In the following sections all these relations will be treated as fundamental, 
and only their consequences will be considered. 

3. Nuclear Physics 

In this section we shall indicate the way in which the fundamental 
parameters determine: (1) which subatomic particles are stable under 
normal conditions; (2) how much energy can be released by thermonuclear 
processes; and (3) the value of the threshold collision energy above which ther- 
monuclear interactions can take place. 

The first of these questions can be answered fro a knowledge of the 
energies of the relevant particles and of the microphysical conservation 
law which affect them, since under low energy conditions (i.e. low 
temperatures, pressure, and velocities) a particle will be stable if and 
only if there is no particle of lower energy to which it can decay without 
violation of the conservation laws. Many quantities are known which are 
approximately conserved in the sense that they can only be changed by the 
action of weak forces and which therefore give rise to particles which are 
stable in relation to the moderately long timescales determined by these 
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forces. However the only quantities which are at present beheved to be 
absolutely conserved are electric charge number, baryon number, and 
lepton number (there may be two distinct and separately conserved lepton 
numbers associated with electrons and muons respectively^^'^^) apart from the 
dynamical quantities - energy, momentum, and angular momentum, whose 
conservation follows from Lorentz invariance. 

The lowest lepton state is the neutrino, which is therefore stable. 
(If muon and electron neutrinos exist as distinct particles, then they are 
both stable by their separate lepton number conservation laws. However 
the electron is also a stable lepton since it is the lightest electrically 
charged particle. The other known lepton, the muon, decays into an electron, 
a neutrino, and an antineutrino; moreover no stable bound states of these 
leptons can be expected, since the strong forces do not affect them, the 
electromagnetic forces between them are always repulsive (both electrons 
and muons have negative charge) and the weak forces, in addition to being 
literally weak, have very short range, and so could only hope to bind 
extremely massive particles, which the lepton, as their name suggests, are 
not; 

None of the above conclusions depends critically on the values of 
any of the microphysical parameters, but this is no longer true when we move 
on to consider baryons. The lowest baryon multiplet is that of the 
nucleons, of which there are two, one charged, i.e. the proton, and one 
uncharged, i.e. the neutron. The mass difference between the nucleons and 
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the next lowest baryon multiplet (the A particle) is very large, 

and therefore the nucleons are the only baryons whose stability need be 

considered except under very high energy conditions. However the 

stability of the nucleons themselves must be examined more carefully 

since their mass splitting, An, (given by (4) ) is relatively 

small, and since there exist reactions by which one could change into 

the other without violation of charge or lepton conservation, i.e. 

n p+ e'+u (/3-decay) and p n+ e'^+u (inverse 

/9-decay). The importance of (8) is now immediately apparent. 

The positive sign of An guarantees the stability of the proton, and 

the fact that it is somewhat greater than me, means that the neutron 

has sufficient energy to be unstable to /9-decay (the neutrino rest mass 

being negligible, probably zero). We see that because An is greater 

than rric, by (8), a neutron will normally be unstable because it can 

decay into a proton. If it had happened that we had me > An > — mg 

then both protons and neutrons would have been stable, while if An 

had been less than —me, then the proton would have been unstable, 

with the result that hydrogen would have been unknown in chemistry. We 

can rewrite (8) in the form 

An - me ^ me 

Because the right hand side is smaller than most masses that occur in 
nuclear physics, only a small circumstantial increase in the energy of a 
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proton, resulting from the effects of its environment, will be sufficient 
to make neutron decay energetically unfavourable, so that a neutron may 
then be stable; the fact that the right hand side is precisely me 
means that if protons and electrons exist together then the same energy 
that is required to make the electrons relativistic will also be just 
sufficient to make it energetically favourable for them to combine with 
the protons to form neutrons. We shall see one of the consequences of 
this in the next section. 

Unlike the leptons, the baryons have many stable bound states, i.e. 
the atomic nuclei. It is now generally accepted that the dominant force 
between nucleons results from the exchange of pions with pseudo-scalar 
coupling governed by the strong coupling constant gs given by (3) 
(heavier mesons, e.g. the K-particles, no doubt contribute, but their 
effect can be expected to be considerably smaller). By the equivalence 
theorem of relativistic field theory^^, such a coupling is approximately 
similar to a pseudo vector coupling with the reduced coupling constant 
i(m7r/mN)5's • Since we only wish to consider orders of magnitude 
we can make a further simplification and treat the nuclear force according 
to the original Yukawa scalar coupling theory according to which the 
potential energy of interaction between a pair of nucleons separated by 
a distance r is gy exp(— r) where gy is the 
Yukawa coupling constant which is of the same order of magnitude as the 
pseudo vector coupling constant. Empirically it is found^^ that the effect 
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of the approximations which have been made is best allowed for by 
dropping the factor 1/2 which leaves the formula 

9y ~ 14 

which gives gy ~ 4/7 or gy ~ 3/10. 

One can treat the two nucleon bound state by an analogue of the 
standard hydrogen atom calculation, provided that relativity can be 
neglected and provided that the internucleon separation is small compared 
with so that the exponential factor in the potential 

can be ignored. The total binding energy in the ground state is then given 
by J gy-fT^'H , from which the binding energy fraction, 
is obtained as ^ 

£n~-^y (15) 

The mean kinetic energy contribution is the same as the binding energy and 

therefore also gives the kinetic energy per unit mass. It can 

be seen from the numerical values that £n is much less than unity, 

and therefore the assumption that relativity can be neglected is 

justified. The mean distance from each nucleon to the centre of mass in 

the ground state will be given by the corresponding Bohr radius, and so the 

mean internucleon distance, will be twice this, i.e. 

Tm ~ 2(7y^'^n ^ • Therefore the mean value of the factor 

in the exponential is rmfn^, ^ 2gY'^m.„ /m^ 2g^'^m^/'mT, by (14), 
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and as a consequence of (7) we see that this is of order unity. 
This means that it is not a very good approximation to neglect the 
exponential factor, with the result that (15) rather overestimates 
the order of magnitude of the binding energy. As soon as higher states 
than the ground state are considered, the exponential factor becomes 
really important, and in fact the nuclear force is so much reduced at the 
higher radii involved that higher bound states cannot exist at all. Thus 
we now see an important qualitative consequence of the coincidence (7): 
if the right hand side had been somewhat smaller many bound states would 
have been possible as well as the ground state, while if it had been 
somewhat larger no bound states would have been possible at all; as it is 
there is just one bound state. 

We can go on to treat the many nucleon problem by the Fermi gas 
model, making the assumption, suggested by the two nucleon calculation, 
that the mean internucleon separation is about . (The 
correctness of this assumption i.e; the fact that the mean separation 
does not change very much with the nucleon number, results from complicated 
saturation effects whose details are still not fully understood). This 
implies that the mean momentum of the nucleons is about rriT^, and, in 
small and medium sized nuclei where electromagnetic effects are unimportant, 
that the potential energy of each one is roughly . Thus the 

mean kinetic energy per nucleon is about ^ m^m^^ (in the 
non-relativistic approximation which is justifiable as before) and the 



19 



mean potential energy per nucleon (dividing by 2 so as not to count each 



which, again by (7) and (14), is the same as the binding energy fraction 
given by (15). Of course subtraction of rough estimates which he 
close together is a dangerous procedure, and ought really to be used only to 
obtain upper limits on absolute magnitudes. Fortunately, however, more 
accurate calculations (supplemented by experiment, because the situation 
is too complicated for pure theory) show that the subtraction is justified 
in this case^^. (in fact both the potential and kinetic contributions have 
been somewhat underestimated, so their difference is nearly correct). 

Thus (15) gives a good estimate of the binding energy fraction of 
the more tightly bound medium sized nuclei from helium 4 onward, even 
though it overestimates the binding energy in the two nucleon case. 
Using (14) we can rewrite it in the form 



Inserting the numerical values we obtain the familiar result ~ 10~^. 
The fact that this fraction is approximately independent of the number of 
nucleons involved means that nearly all the energy available will 
be released in the first stage (i.e. burning of hydrogen to helium) of 
thermonuclear evolution. Subsequent stages involving the burning of helium 



interaction twice) is rather less than | Qy m.^/m^ . By (7) 

and (14) the former is just half the latter, and soon subtracting we 

find for the total binding energy fraction roughly | gy m.,^/m^ 




(16) 
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to form heavier elements will give relatively negligible energy production. 

In the light and moderately heavy nuclei we have been considering 

so far there are approximately equal numbers of protons and neutrons in 

stable states, because this ratio minimises the effects of the exclusion 

principle and because (as has already been remarked) the mass difference 

An — rric is small compared with the other energies involved. 

Thus we shall have Z ^ ^ Nj^ where Z is the charge number and A^a 

the total nucleon number of the nucleus. 

Since the mean separation of nuclei is about m~^, the nuclear 
radius is of order ^ ^iid so the electrostatic potential 

1 /3 

energy Ec is given in order of magnitude by Ec ~ 2Z^ rriT^ N^^ 
Putting Z ^ ^ Na gives 

E.^^e^m^Nl^' 

from which we see that the electrostatic potential energy fraction is 
of order | e^{mT^/m-^)N^^^ which becomes comparable with £n 
as given by (16) when A'a ~ 60. As A'a increases above 
this value the electrostatic repulsion steadily reduces the binding 
energy fractions and the charge number Z becomes considerably reduced 
below the optimum value | A^a determined by the exclusion principle; 
finally stable nuclei are prevented from existing at all beyond A^a ~ 240. 
Due to the relatively low energies involved, transmutations between heavy 
nuclei are of only minor importance in nature, although 
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decays of heavy nuclei are responsible for heating the earth's interior 
sufficiently to produce a molten core, and to provide the energy source for 
earthquakes and volcanoes; it is also possible that such decays are needed 
to produce some of the heavy trace elements which seem to play an 
important part in biological systems. 

We shall need one more result from nuclear physics, namely the 
threshold temperature for thermonuclear reactions. At low temperatures 
nuclei are prevented from interacting primarily by the electrostatic 
repulsion barrier: two nuclei will be able to react only if they collide 
with sufficient energy to surmount or penetrate it. The electrostatic 
potential energy which has just been calculated can also be used 
as a good estimation of the order of magnitude of the height of the repulsion 
barrier between a pair of nuclei of the corresponding size. It is apparent 
that the barrier is lowest when small nuclei are involved. Such nuclei 
will react most rapidly with each other when their mean kinetic energy is 
comparable with the barrier height, and by (12) the temperature Tn 
at which this occurs is given by 

Tn ^ e^m^ . (18) 

However we should expect the minimum temperature for slow thermonuclear 
burning to be considerable lower than this, since in a Maxwellian 
distribution some of the nuclei will have very much greater energy than 
the mean, and since in addition there is the possibility that a particle may 
penetrate the barrier by quantum mechanical tunneling, even when it has 
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insufficient energy to surmount it. Accurate calculations^^ show that 

slow burning from hydrogen to helium becomes important at a temperature, Th 

which is less than Tn by a numerical factor of order 10~^, i. e. 

we have 

Th ^ IO-^Tn ^ IQ-^e^ mN . (19) 

The calculations are complicated because the process can only proceed 
indirectly via a chain of reactions: either the proton-proton chain 
or the C-N-0 chain. Both chains involve /3-decays(which result from 
weak forces) and consequently the rates are limited even under the most 
favorable circumstances. The density also affects the rates, not only 
by determining the concentrations of the constituents, but also because 
electron screening of the electrostatic barriers can increase the rates 
at high densities. Moreover the C-N-0 chain rate depends on the 
concentration of medium weight elements. However the rates of both 
chains are very temperature sensitive (the C-N-0 chain extremely so) 
and as a result (19) is a good estimate of the order of magnitude 
of the threshold temperature under practically all natural circumstances. 

In a similar way, the temperature Thg for the burning of 
helium to oxygen and carbon is given by 

Thc ~ IO^^Tn ^ lO^^e^ mN , (20) 

while the medium weight elements burn at a temperature of order 10~^Tn 
to form the most tightly bound nuclei with A^a ~ 60, particularly 
ion 56. 
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At temperatures of the order 10~^Tn and below, all the important 
thermonuclear reactions are exothermic? However as the temperature becomes 
comparable with Tn itself, the situation changes radically. Not only 
are there no more exothermic reactions available, after the medium weight 
elements with A'a have been transformed to the most tightly 
bound elements with A^a ~ 60, but in addition endothermic 
reactions reversing the process become possible. The reason for this is 
that, as has been already remarked, the differences in binding energy fractions 
between helium-4 and the medium weight and most tightly bound elements 
are small, in fact by factors of order 10~^ compared with the total binding 
energy fraction e-^. Since the electrostatic repulsion energy 
fraction becomes comparable with e-^ itself for the most tightly bound 
nuclei, it will be comparable with 10~^£:n even for light nuclei. In 
other words Tn (defined as the electrostatic repulsion energy between 
light nuclei) is also the energy required per nucleon for dissociation of 
tightly bound nuclei into helium-4 nuclei. Therefore at temperatures not 
far below Tn, where thermonuclear interactions can proceed most freely 
without electrostatic inhibitions, fusion reactions will cease to be 
energetically favourable, and will be replaced by endothermic dissociation 
of medium nuclei into helium. Because of this drastic change of conditions, 
Tn marks the boundary of the realm of high energy astrophysics. 

We now notice some important consequences of the coincidences (6) and 
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(8), which have the result that two other completely independent effects 
become important at this same temperature. The first of these effects is 
electron positron pair creation which occurs on a large scale at temperature 
Tp which is of course the same as the temperature at which electrons 
become relativistic i.e. we have 

Tp ^ me . (21) 

The second of these effects is the Urea process whereby electrons interact 

with neutrons in the nuclei producing neutrinos, according to the cycle 

e~ +p — > n + z/ followed by n — > p+ e~ + u . During 

the intermediate stage the charge number, Z , of the nucleus involved 

changes bu one from the minimum energy value, which (by a simple order of 

magnitude calculation, using the Fermi gas model) can be expected to reduce 

the total binding energy of the nucleus by an amount of order ^n^^n^a"^ 

in the case of a nucleus where Nji^ is odd, by the sum of two amounts 

of this order where A'a is even and Z is even, and by the difference 

of two amounts of this order when Nji^ is even and Z is odd. In 

order to have a good chance of taking part in such a reaction, an electron 

must have enough kinetic energy not only to make up for this loss of binding 

energy in the nucleus, but also to make up the extra energy. An — me, 

required for the combination with the proton to form the neutron. The binding 

energy requirement is at most of order ej^nif^N^^ and this 

becomes less than An — ttLc for nuclei with A'a greater 

than about 20, so that the latter quantity gives the basic energy requirement 



25 



for the Urea proeess. Therefore the minimum temperature Tu at which 
the Urea process can proceed rapidly should be given by 



actually the Fermi gas model is too crude for satisfactory calculations of 
the energy differences between nuclear states, and in fact the Urea threshold 
depends rather critically on the individual nuclei involved^^; this formula 
is roughly correct for certain particular medium weight nuclei, but if 
these nuclei happened to be absent the minimum temperature would be 
considerably higher. The astrophysical importance of pair creation is first 
that by absorbing kinetic energy it decreases the adiabatic index, 7 , 
of the gas, thus reducing its resistance to compression, and second that 
it, like the Urea process, produces neutrinos, as a result of pair annihilation 
e"*" + e~ z/ + z/ . Neutrino production by either of these processes 
is an endothermic reaction which is normally irreversible, since, with their 
very low interaction cross sections, neutrinos will usually escape easily 
from any limited volume of gas. Pair creation is rather less temperature 
sensitive than the Urea process, and can have a significant effect considerably 
below the characteristic temperature i.e. at temperature of order 10~^Tp. 

As a consequence of the microphysical coincidence (6) we can now 
derive from (18) and (21) the astrophysical coincidence 



Tu An - rric 




(23) 



and similarly as a consequence of the form (13) of (8) we can derive 
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Tu^Tn. (24) 

This double coincidence means that the endothermic processes resulting 
from thermonuclear dissociation, pair creation, and the Urea process all set 
in within the range 10~^Tn to , a fact which considerably 
complicates the analysis of high energy phenomena, such as the catastrophic 
collapses or supernova explosions which result in a star that reaches these 
temperatures. It is under these conditions that the heavy elements, with A'a 
greater than 60 (which cannot be produced exothermically) are believed to 
have been created. 

4. The Equation of State of Cold Matter 

At sufficiently low pressure, cold matter exists in a solid or liquid 
state, in which the density is determined essentially by a balance of 
electrostatic attractive forces (between positively charged nuclei and 
negatively charged electrons) and repulsive forces of quantum mechanical 
origin. It is sometimes stated that the balance is achieved by electrostatic 
repulsion between electrons, but this is rather misleading, since there are 
equal numbers of positive and negative charges, and so the particles will 
be able to arrange themselves roughly in a lattice in such a way that 
between nearest neighbours attractions are dominant; in classical 
theory collapse would inevitably follow, but in quantum physics it is 
prevented by the exclusion principle, since electrons are fermions. The 
net effect is most easily calculated when only two particles are 
present, one of each sign, i.e. in the case of a single hydrogen atom. 
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The solution of this famous problem gives the mean separation of the proton 
and the electron in the most tightly bound state as the Bohr radius, e~^m~^, 
the corresponding binding energy being | e^me .. Heavier elements 
form atoms whose radii are still given in order of magnitude by the Bohr 
radius '(since although there bare many electrons present to maintain charge 
neutrality, the inner ones are correspondingly much more tightly bound due 
to the stronger attraction of the nucleus), and the binding energy of the 
outer electrons will still be given in order of magnitude by | rUe 
(for the inner ones it may be considerably larger). When many nuclei are 
present the situation is not very different, except in degree of complexity; 
they will associate in molecules or crystals whose density is only a little 
less than that of the separate atoms, i.e. the mean separation between nuclei 
will still be a few times ■ Thus a typical nucleon number 

density riN , in a solid or liquid at low pressure will be given by 

riN ^ X 10"^ Na m ^ [25) 

where A'a is the mean nucleon number. (The precise value 1/7 

in the numerical factor has been chosen so that the result is almost exactly 

correct in the case of water (with A^a = 6) and iron (with A'a = 56). 

the density would be rather higher if it were not for the highly directional 

nature of chemical bonding effects, which tend to produce many empty 

spaces in the structures (as a consequence of which water can absorb large 

quantities of suitable solutes without much change in volume). As a result 
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of the association the binding energy of the most loosely bound electrons 
will be somewhat increased, typically by a factor of order 10^^ . It is 
these marginal changes in binding energy that are important in chemical 
rearrangements^^. Thus we may estimate a typical order of magnitude of the 
chemical binding energy fraction as 

ec-^x 10-\mjNAm^) (26) 

We may compare this with the nuclear binding energy fraction (16): thus 

we have, in the case of the lighter elements, Ec/^n ~ \{e/gsyi^N/^n)^{^c/'m'Tr), 

or numerically ^c/^n ~ 10"^ , the smallness of this being mostly 

due to the factor {e/ g^Y since the mass ratio factors cancel out. 

(This shows incidentally why a rather inefficient hydrogen bomb weighing 

a few tons can be compared with a T.N.T. bomb weighing several millions of 

tons.) In absolute terms we see that Eq varies from about 10~^ for 

light elements to 10^^^ for heavy elements. 

Using the energy values that have been calculated, we can evaluate 
some important temperature limits for use later on. First Tj , the 
temperature at which light elements must be fully ionised, must be given by 

TiRi^e^me. (27) 

Actually hydrogen will be partially ionised well below this temperature since 
in a Maxwellian distribution some atoms will have much more than the mean 
energy and so will be able to cause ionisation by collisions. The temperature 
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at which hydrogen becomes sufficiently ionised to be highly opaque is of 
order 10~^Ti , and significant opacity first appears even below 
this. Second Tc , the upper limit of chemical temperatures, above which 
all molecules will be fully dissociated, must be given by 

Tc ^ ^ X 10^^ me (28) 

which is the same as the temperature of partial ionisation of hydrogen. 
Again this is an extreme upper limit. Most crystals and large molecules will 
break up at temperatures of the order of 10~^Tc or lower. Finally we may 
estimate Tb , the temperature at which biological processes take place. 
The basic structures of important biological molecules - proteins and nucleic 
acids - depend on strong chemical bonds, but the 3-dimensional arrangements 
which they take up (which are vital for their specific functions) are 
determined by weak hydrogen bonds between different parts of the same or 
neighbouring molecules. Hydrogen bonds^^'^^ result from small electrostatic 
perturbations of the basic chemical structures, and as a result they are 
weaker than ordinary chemical bonds by a factor of order 1.20. At 
temperatures of the same magnitude as the hydrogen bond energy, proteins 
become denatured (i.e. they lose their characteristic 3-dimensional 
structure) and water (whose molecules are also held in contact with each 
other by hydrogen bonds) boils, so that biological systems are destroyed^^. 
However at much lower temperatures the hydrogen bonds become in effect 
completely rigid, so that the continual rearrangements that are the essential 
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characteristic of life processes are unable to take place. Therefore 
biological systems can function only in a unique temperature range, which 
is less than, but not much less than, the typical hydrogen bond energy, i.e. 
in the neighbourhood of a temperature Tb given by 

Tb f« lO-^e^me. (29) 

This is the fundamental formula for what is commonly called room temperature. 

All the temperatures that have just been calculated give limits on the 
conditions that can be considered to be cold in various contexts. For our 
purposes in this section "cols" means sufficiently cold for matter to be 
in a solid or liquid state in order that the formula (25) can be applied. 
Most substances satisfy this condition at a temperature around Tb , 
but some exceptionally tightly bound and symmetric molecules remain gaseous 
well below this temperature - most notably the helium atom which is spherically 
symmetric in its electronic structure when isolated (and so almost unaffected 
by electric perturbation fields), and which remains gaseous^^ down to temper- 
atures not much greater than 10~^Tb. Despite of this, the formula (25) 
can often be used at temperatures comparable with Tq , for example when 
heavy elements are predominant, or when there are moderate pressures present 
(but not sufficient to cause a high degree of compression of the solid or 
liquid) - both of which conditions are satisfied in the centre of the earth. 
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We can now move on to consider how the situation is modified when 
high pressures are apphed. The expression (25) gives a vahd estimate 
of the density so long as the pressure forces are small compared with the 
electrostatic attractions between the nuclei and the electrons. However 
for larger pressures we can ignore the electrostatic effects and balance 
the repulsive effects due to the Pauli principle directly against the pressure, 
or in other words, we can treat the cold matter as a degenerate Fermi gas, 
following Chandrasekhar (1939)^. Electrons will no longer be associated with 
particular nuclei (an effect often described as pressure ionisation) so that 
the situation is much simpler than at low pressure. If there are n 
fermions of a particular kind per unit volume, then their mean momentum, p, 
will be given approximately by p n^/^ . If each fermion has mass m 
then the magnitude P of their contribution to the pressure is given by 
P ^ np or P^np'^/m according to whether the particles are 
relativistic or not; i.e. according to which formula gives the smaller 
contribution. Thus the degenerate Fermi gas pressure is P ~ n^^^ or 
P m~^n^^^ , whichever is smaller. 

At moderately low pressures the non-relativistic formula will be the 
relevant one, and since the mass occurs in the denominator the lightest 
particles, i.e. the electrons, will give the dominant contribution. The 
electron number density Uc can be estimated as rig ~ {Z / Np^n^. 
We have remarked that the factor Z/N^ varies from 1 for hydrogen to about 
I for heavier elements, so we can ignore it in a rough calculation and use 
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He ^ n-^- Thus we obtain the formula 

T-> -1 5/3 

P n^' . 

It will become valid when n-^ becomes larger than the ordinary solid state 
density. Using (25) we see that the critical pressure Pc at which this 
crushing sets in is given roughly by 

Pe-lO-'iV^'/'m^V". (31) 

Below this value the density is approximately independent of the pressure. 
Actually this pressure can be expected to be an overestimate since there must 
be a transition region where both electrostatic attraction and pressure are 
effective, (the transition region will be longest for heavy elements). Exact calcu- 
lations for iron have been made by Feynman, Metropolis, and Teller^^. (See also 
Knopoff)^^ 

The formula (30) remains valid until the density 

nN ~ (32) 
is reaches, corresponding to the pressure 

P ^ mt (33) 
at which stage one might expect at first sight that the relativistic formula 

P ^ n^^^ (34) 

would take its place. Actually however, complications arise due to the 
coincidence (8) (or more transparently, due to the alternative form (13) ) 
by which the point at which the electrons become relativistic is also the point 
at which they have sufficient energy to undergo inverse /3-decay. Thus at 
the density and pressure given by (32) and (33) the electrons will begin 



33 



to combine with the protons in the nuclei to form neutrons. Since the extra 
neutrons will tend to be lost from the nuclei, Wheeler^'' has described this 
situation as the regime of neutron drip. The net effect of this is that 
the the pressure will drop substantially below that given by (34). Eventually, 
when the pressure is raised sufficiently, the matter will become effectively 
a neutron gas (the electrons and remaining nuclei giving an insignificant 
effect) so that the formula 

P ^ m^^ n^/' (35) 
can be used. This will probably remain valid until the transition density 

rzN ~ (36) 

corresponding to the pressure 

P^m^ (37) 

at which the neutrons themselves become relativistic, when one might expect 

a return to the formula (34). There is considerable uncertainty at this stage 

however due to the fact that in between the density (33) at 

which electrons become relativistic, and the density (36) at which the 

neutrons become relativistic, is the nuclear density, which by the results 

of section 3 is 

nN ~ . (38) 

Present day understanding of nuclear forces is insufficient for making 
reliable calculations much beyond this density, but it is thought that 
repulsive effects may become important, in which case (43) will be an 
underestimate of the pressure. However very divergent opinions about this 
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have been held by different authors -Cameron (1959)^^''^^ favoured an extreme 
hard core repulsion between nuclei (in his original paper he overlooked 
the fact'^'^ that causality in special relativity requires that the hardness 
must be limited), while Ambartsumian and Saakyan (1960)'^^ favoured a 
very soft equation of state postulating that the nucleons can evade the 
effects of the exclusion principle by changing into heavier strange particles 
such as the A hyperon), and Harrison and Wheeler (1964)^° have favoured 
a middle way by which (25) is extrapolated without modification. This 
situation has been reviewed by Salpeter (1965)^^ and also by Zel'dovich 
and Novikov (1965)^^. Fortunately, in considering the equilibrium of 
astrophysical bodies, it turns out to be unnecessary to go much beyond the 
density given by (36), as will appear in the next section. 

5. The Equilibrium of a Cold Spherical Body 

We shall first consider the general problem of a spherical body 
containing nucleons whose gravitational self attraction is balanced by a 
pressure gradient. Provided that the body has a simple density profile, very 
simple considerations enable one to make good order of magnitude calculations 
of the conditions for equilibrium. By a simple density profile is meant one 
where the density decreases to a small fraction of its central value at a 
more or less well defined radius R , such that only a correspondingly 
small fraction of the total mass lies further out than R . An example 
of a more complicated density profile would be one where a dense core is 
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surrounded by an envelope, containing the greater part of the mass, but 
with a density nowhere exceeding a very small fraction of the typical core 
density, as is the case for example in highly evolved stars which have 
reached the red giant stage of evolution^^. However accurate calculations 
show that the assumption of a simple density profile is satisfactory in all 
the cases that will be considered in this and the subsequent sections, and 
it will be taken for granted from now on. Thus if tt-n is the mean nucleon 
number density, then the radius R at which the density drops to a 
negligible fraction of the central density will be given by 

R ^ n^^^^N^/^ . (39) 

The mass, M, is of course given by 

M^Nm^. (40) 

We shall use Newtonian gravitational theory, but we shall need to check that 
this is adequate in each case we consider by verifying that General Relativity 
effects are small. For this purpose we shall need the formula 

M / R ^ n^^^ N"^^^ . (41) 

So long as M/R is small. General Relativity effects will indeed give 
only minor corrections (although in some critical situations even these minor 
corrections will be found to be important). However when we have M/R r^i 1 , 
the body will be lying close to its Schwarzschild radius'^", and Newtonian theory 
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will begin to fail completely. By (41) there will be no danger of this 
occurring so long as the density is less than the critical value given in 
terms of the nucleon number by 

Un ^ m^^ N-^ (42) 

but if the density becomes larger than this, General Relativity effects will 
become predominant. In principle the body might then balance in unstable 
equilibrium'^'^, a little outside its Schwarzschild radius, but in practice it must 
either expand back to lower densities or else undergo irreversible gravitational 
collapse within its Schwarzschild radius. 

In order for a body to be in equilibrium, the gravitational force density 
must balance the pressure gradient. The typical mass density is mN,?T-N 
and so a typical gravitational force density is of order mN M/ E? . 
A typical value for the pressure gradient will be of order P/ R where P 
is the value of the pressure at the center. Thus by (39) and (40) the 
equilibrium condition must take the form 

which is just a simple version of the virial theorem. 

We can now use this formula in conjunction with the equation of state 
given in the previous section to determine the mean equilibrium density of 
a cold spherical body as a function of its nucleon number. 

Provided that the central pressure P is less than the transition 
value given by (31), a cold body will be supported by ordinary solid or 
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liquid state forces at the density given by (25). This condition will 
hold so long as is less than the critical crushing value, Nc say, 
which is obtained by substituting (25) and (31) into (43). Thus 
we have 

N^^^Nl^'e'N^, (44) 

where by (2) the Landau number, A^'l, is defined by A^'l = m^^ . 
Similarly in terms of the Landau mass Ml = , we obtain the 
corresponding critical mass as 

Beyond the crushing limit Nc, serious compression sets in and the 
degenerate gas formula (30) must be used. For a body such as the earth, 
which is made largely of silicon and iron with a relatively small contribution 
from the lighter elements^^, we may use A^a ~ 30 as a rough average. 
The earth (whose nucleon number is ^ 2 x IO^^A'^l ) is in fact 
below the crushing limit but not very far - actually within a factor around 
10^^ . Jupiter, which is 300 times more massive, lies actually at the 
crushing limit, which is rather smaller in this case, since being made 
largely of hydrogen, with a small contribution from other light elements 
such as carbon; Jupiter has A^a ~ 4 . Since the limit is not sharply 
defined, this means that Jupiter is in fact compressed considerably 
beyond normal densities. 

Above the crushing limit the density is given in terms of the nucleon 
number ( (30) and (43) ) as 



nN ^ m^iN/Ni^)^ 



(45) 
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A cold crushed body in this range is in fact a black dwarf star. However 

for the higher masses of the range where the degeneracy pressure is fairly 

high, such a star could have a correspondingly high temperature without much 

effect on the equation of state, and so (45) can also be used to 

describe white dwarf stars^^. This formula will be valid until the electrons 

become relativistic, at the density given by (32) which we see occurs 

at the critical nucleon number 

N ^ Nj^. (46) 

This critical point (originally predicted by Landau) was first studied 
in detail by Chandrasekhar7'^ 

If it were correct to use the relativistic formula (34) at these 
densities we would find that beyond this stage the density would be 
indeterminate: the matter would exist in a state of metastable equilibrium 
with remaining at the value TVl. However we have seen that at these 
higher densities the pressure actually drops below that given by (34) as a 
result of the neutron drip effect (which occurs at this point in consequence 
of the coincidence (8) ), and the effect of this is that the equilibrium 
nucleon number actually decreases as the density is raised. An equilibrium 
of this sort is obviously unstable. 

The nucleon number rises again when the formula (35) becomes valid, 
so that by (43) the density will be given by 

riN ^ {N/N^f . 47 
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A body in this range is a neutron star. It is clear that it can have a very 
high temperature, of order rric or more, without much affecting the equation 
of state, i.e. while remaining effectively cold as far as the equilibrium equation 
(47) is concerned. Such a body can have a correspondingly high surface 
temperature, and is a potential X-ray source. '^^'^^ We see that a neutron star 
is denser than a cold dwarf star of the same nucleon number by a factor 
[m-^/me)^ which numerically is of order 10^°. 

If nuclear forces were unimportant, the formula (47) would 
remain valid until the density reached the critical value given by (36), 
when the neutrons would become relativistic, after which the formula (34) 
would become valid again and as before there would be a state of metastable 
equilibrium with N ^ Li^ . We should then deduce that no cold equilibrium 
state is possible at all with much larger than Al '■ such a body 
would be forced either to explode or to collapse. 

Despite of the fact that nuclear forces probably do have an important 
effect on the equation of state, this last conclusion is almost certainly 
correct. The reason for this becomes apparent as soon as we check, as we 
must, that Newtonian theory is valid for the calculations we have made so far. 
We easily verify from the formulae that have been obtained that, at all stages 
short of the neutron star limit, the equilibrium lies well below the General 
Relativity limit (42). However at the stage when degenerate neutrons 
become relativistic we have riN ~ (by (36) ) and A^ ^ Al = , 
and on substituting into (42) we see that the neutron star limit lies 
just at the point where General Relativity first becomes important. This was 
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first realised by Oppenheimer and Volkoff, wlio made tlie original study of tliis 

limit (It is worth emphasizing that this effect holds automatically - it 

is not a coincidence of a coincidence such as (6), (7), or (8).) 

As a result of this, it does not matter very much that the equation of state 

is uncertain at densities higher than , since a body will collapse 

within its Schwarzschild radius before such densities are reached. 

(Due to the effects of space-time curvature it is possible in principle for 

a body to exist in equilibrium near to its Schwarzschild radius at much 

higher densities^^. However such equilibrium is necessarily unstable, and 

so can never be maintained in practice.) Unfortunately the equation of state 

is uncertain rather below this density (since the nuclear density riN ~ 

is lower than this by a factor of order 300). The effects of this are not 

important for order of magnitude calculations, but exact calculations are 

upset sufficiently for it to be still a matter of controversy whether the 

exact value of the Oppenheimer - Volkhoff nucleon number is slightly more or 

slightly less than the Chandrasekhar limiting value. This is a pity, because 

neutron stars are likely to be extremely rare if the latter is true, unless 

it is possible for one to be formed from a star originally larger than the 

Chandrasekhar limit which collapses beyond it, and after doing so manages to 

explode as a supernova leaving its core behind as a neutron star'*'*'^^'^^. 

The timing is critical, because if the explosion occurs too late, the core will 

already have collapsed within its Schwarzschild radius. 
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6. The Equilibrium of a Hot Star 

By a hot star we shall mean a spherical body where the central 
temperature is so high in relation to the density that the effects of degeneracy 
and of valence, electrostatic, and nuclear force can be ignored. In such a 
body the equation of state will be that of a non-relativistic perfect gas in 
equilibrium with electromagnetic radiation, i.e.® 

P = nT+ —T^ (48) 
45 ^ ' 

where n is the number density of gas particles. 

We shall work separately with two distinct cases which arise when 
either one of the gas or the radiation contributions dominates the other. In 
order to distinguish the two we shall make use of the coefficient i] defined as 
the ratio of the gas contribution to the radiation contribution. (This will 
be rather more convenient than the conventional coefficient f3 defined as 
the ratio of the gas contribution to the total pressure; they are related by 
7] = — P) .) Thus we have 

V 

and we have the pressure laws 

P 

when ?7 ^ 1 , and 

P 

when ?7 <^ 1 . 



nY (50) 

(51) 
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We shall only use these formulae for stars with mean temperatures in 
the range between the temperature of total ionisation, Tj ^ ^e^ rrie , 
(by ((27) ) and the temperature at which pair creation becomes important, 
Tp ^ TTie (by (31) ). Within this range we shall have n ^ . 

(When T is greater than Tp , n will be considerably greater than 
riN , due to the creation of positron electron pairs, and, at even 
higher temperatures, of mesons etc; (and in any case the electrons will 
cease to satisfy the non relativistic assumption) but no star can exist 
for long at these temperatures without a catastrophe in the form of an 
explosion or coUapse^'^'^^. When T is less than Tj ^ ^e'^ , the gas 
will be at most partially ionised, and n may drop considerably below riN 
depending on the atomic weight of the atoms present. However a star 
cannot remain for very long at such a low mean temperature, since it will 
have a relatively low opacity also, and will therefore lose energy rapidly 
by radiation; if it is nearly degenerate it will simply cool and turn into 
a black dwarf as a result; otherwise it will contract until either degeneracy 
sets in or until the temperature is raised above Tj.) 

It will sometimes be desirable to use a more accurate formula than n/riN ^ 1 
when the factor n/riN appears as a high power, and in such circumstances 
we shall use n/riN ~ 2 , which will be correct when the star under 
consideration consists mainly of hydrogen, as it will during the greater 
part of its evolutionary lifetime. 

When 7] ^ 1 we have, by (43) (the virial theorem) and (50), 
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the equilibrium condition 

riN ~ 8(iV/A^L)"^T^ (52) 

and hence, by (49), 

r^^W{N/N^)-\ (53) 

Similarly when -C 1 we have by (43) and (51) the equilibrium 
condition 

n^^^-{N/N^y'''T' (54) 

and hence, by (49) 

T^^?,{N/N^)-^'\ (55) 

We see that the ratio rj depends only on the mass, and that the transition 
region r/ ~ 1 occurs at the critical nucleon number A^^ given by 

AT, ^ lOiVL . (56) 

Thus we have found a second significance for the Landau number: for stars 
up to ten times this size gas pressure will be dominant and the formulae (52) 
and (53) will be valid; while for stars larger than this, radiation pressure 
will be dominant and the formulae (54) and (55) will apply.*'^^ 

At temperatures greater than Ti, a star will normally lose energy 
steadily by radiation (despite the relatively high opacity). As a result 
it will in general contract, continually raising its temperature in accordance 
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with the virial relation, T oc n^^ , which holds whether gas or 
radiation pressure is dominant, since it follows from both (52) 
and (54) 

In the case of a star below the Chandrasekhar limit, the 
contraction may be halted if it becomes degenerate, in which case, in this 
temperature range, it will become a white dwarf. (It will be white because 
as a result of the extremely high conductivity in the degenerate region 
where the material is effectively metallic, the energy losses will be 
limited only by the opacity in the surface layers, so the surface temperature 
will be much higher than in a normal star with the same mass 
and temperature.) Radial contraction will then be replaced by a steady 
decrease in luminosity as the star cools toward the final dead black dwarf 
state^^. 

However contraction will also be halted if, instead of becoming 
degenerate, the star reaches a sufficiently high temperature for hydrogen 
burning to take place a a sufficient rate to replace the energy lost by 
radiation. In this situation, the main sequence stage, it will be able to 
remain with a more or less constant radius and luminosity for the greater 
part of its evolutionary life. When the hydrogen in the central regions 
is exhausted, contraction will take place again, leading to higher temperature 
in the range from Th to Tn where burning of heavier elements takes 
place, but no comparably long steady state will result because much less 
energy is available from these reactions. Beyond this stage, if the body 
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is too massive to have a cold equilibrium state, it will have to explode 
or contract. We have seen by equations (23) and (19) that as a 
result of the coincidence (6) the temperature Tn coincides with the 
upper limit Tp of the range in which the above formulae are valid, 
while the temperature Th lies near the middle of the range from Ti to Tp 

The lower limit on the size of a range of hot stars defined with a fixed 
temperature T occurs where the locus in the n^, N plane intersects 
the degenerate electron region: thus by (43) and (52) this lower size 
limit is given by 

iV/iVL ^ {T/m,f/^ . (57) 

It happens (and no coincidence is involved) that at the upper limit Tp ^ me 
of the temperature range under consideration, this intersection coincides 
with the upper limit of the degenerate electron region, that is the 
Chandrasekhar limit point with ^ Ni^ = m^ and tzn ~ . 
Similarly at the lower limit Tj | me of the temperature range 
under consideration, this intersection coincides with the lower limit of 
the degenerate electron region, that is the point at which serious 
crushing sets in with N Nc ^ e^m^^ (in the case of hydrogen with 
unit atomic weight) and tt-n ~ m^ . However of more importance 
than these is the lower limit, A'h say, of the size of a hydrogen 
burning star, with temperature Th ~ 10~^e^m^ by (19), which 
by (57) is 

Nn^lO~'(^] iVL. (58) 
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Thus as a result of the numerical coincidence (6) we obtain finally 
A'^H ~ lO^^A'^L , a value that is confirmed by the more accurate calculations 
of Kumar (1963).^^ This is the absolute lower bound for steadily 
burning stars. It differs from the Landau number only by an arithmetical 
factor, which, it will be remembered, arose from the properties of the 
Maxwell distribution. A star that is below the Chandrasekhar limit by 
more than this factor will reach equilibrium as a degenerate body before 
it is hot enough for thermonuclear burning. 

The upper limit on the range of hot stars with a fixed 
temperature T occurs at the General Relativistic limit: thus by (42) 
and (45) this upper limit is given by 

N/N^ ^ {T/m^f . (59) 

For very large hydrogen burning stars the temperature must be rather larger 
than the minimum value, Th ~ 10~^Tn , given by (19) since the 
required rate of energy output becomes so large that despite the high 
temperature sensitivity of the thermonuclear reactions a noticeable temperature 
increase is required to maintain it. Thus in such stars the value 5 x 10~'^Tn 
(at which the C.N.O. cycle completely dominates the proton-proton 
cycle) is a rather better estimate of the central temperature. The 
difference becomes significant in the formula (59) where the square of 
the temperature is involved. Inserting this value we see that the maximum 
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conceivable nucleon number, A'^q , for a supermassive hydrogen burning 
star is given by 



which, using the numerical values, gives A'^q ^ IO^^A^l • However in 
practice it is difficult for a star to exist with a size anywhere near as 
large as this, as result of the tendency to instability which arises 
whenever radiation pressure is dominant. The reason for this is that 
under adiabatic compression radiation pressure varies with volume according 
to an inverse 4/3 power law, and by equation (43) this leads to a 
metastable equilibrium which can easily be converted to instability. This 
effect becomes important for stars a few times larger than the transition 
size where gas and radiation pressure are comparable, and therefore we 
estimate the upper limit to the size of a stable star as 



Careful calculations by Schwarzschild and harm (1958)^"^ indicate that 
stars will in fact be pulsationally unstable when N is larger than 
about 40 A'l (i.e. 65 Nq) so that this estimate is fairly generous. It 
appears that if a star larger than this starts to condense out of a gas 
cloud, the radiation pressure builds up to such an extent that the outer 
parts are blown away, and only a central core below this size remains. For 
stars with larger than about 10^ A'l this tendency to instability 




(60) 



A^j ^ IO^ATl . 



(61) 
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is reinforced by the marginal effects of General Relativity, which is thus 
able to have an important effect long before it would have been able to in 
a less precariously balanced situation. This is one of the reasons why it 
is difficult to create a simple theory which explains db QUcLScir clS Sb 
supermassive star. (The tendencies to instability might however be 
counteracted by rotational^^ and magnetic effects, so it remains conceivable 
that a slightly more complicates theory of this sort could succeed. However 
it should be remembered that it is by no means certain that quasars are 
supermassive objects at all: the local hypothesis remains very much alive). 

This completes our demonstration of why stellar nucleon numbers must 
lie close to the Landau number. The upper limit A'^j ^ A'^l results 
from the fact that the Landau number determines the boundary between gas 
and radiation pressure, while the lower limit A'h ~ 10~^A'l , 
(which could not conceivably lie above the Chandrasekhar limit) lies close 
below A'l as a consequence of the coincidence (6). The sun lies 
near the middle of the allowed range with A^0 ^ Ni^ . 

The results of the last two sections have been collected together in 
a logarithmic diagram of the n^, N plane. No attempt has been made 
to make the diagram more accurate or realistic by smoothing out corners 
where straight lines meet: it has been deliberately left in its crude form 
in order that its primary structure should not be obscured. 

In the course of its evolution a spherical body will normally move 
horizontally from left to right across the diagram as it loses energy 



49 



Legend for Figure 

Total nucleon number, A^, is plotted logarithmically against 
mean equilibrium number density, riN- 

Loci corresponding to the main equations of state of sections 5 and 6 
are plotted as follows: 

Cold bodies (broken line in the unstable neutron 

drip region). 

— Hydrogen burning stars, T Th, (broken line in 
unstable region). 

— - — - — - — - General Relativity limit. 

Limit of complete ionisation. 

(The position of red giants and of stars in the last stages 
before a supernova explosion have been included, although they are not 
susceptible to the simple analysis given here. The indicated values of 
riN are intended to suggest roughly the densities in the most central 
core, and bear no relation to the outer densities. In the same spirit 
the limiting region T ^ Tp at which pair creation and endothermic 
processes become dominant has been indicated thus: oooooooo.) 

The shaded region indicates the boundary to the right of 
which no stable spherical body can exist. 

The wavy line separates the region where radiation pressure 
is dominant (above) from the region where the gas pressure is dominant 
(below) . 
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by radiation, with a long almost stationary period before it passes through the 
hydrogen burning zone. Its motion may have an upward component if it 
accretes matter, or a downward component if it ejects matter (as in planetary 
nebulae). It may move downward discontinuously if there is explosive 
ejection of matter as in a supernova. Ultimately the body will come to 
rest if it reaches a cold equilibrium position, or undergo gravitational 
collapse if it reaches the General Relativity limit. It can be seen that 
most, and perhaps all, of the neutron star range is inaccessible except to 
a body that loses matter at some stage. 

7. STELLAR LUMINOSITIES 

The luminosity of a star is governed by the rate at which energy leaks 
out from the interior regions, and this in turn is governed by the opacity 
of the matter in the intermediate regions if the energy transport is primarily 
by radiation, and otherwise by the rate of convection. It is sometimes 
stated that the rate of energy loss in a main sequence star is governed by 
the rate of thermonuclear energy generation, but that is a physically 
misleading way of thinking about the problem. The rate of thermonuclear 
energy generation is, as has previously been remarked, an extremely sensitive 
function of temperature, so that our rough order of magnitude estimates of 
the temperature would be quite inadequate for calculating it directly. For 
the same reason it is unnecessary to do so. What happens in fact is 
that the rate of energy generation automatically adjusts itself to balance 
the losses by radiation and convection: if the rate is a little too high an 
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energy surplus will build up in the star causing it to expand, thereby 
reducing the central temperature until the rate of generation has 
dropped to the rate of escape; similarly if the rate of generation 
is too low the star will contract until it has been raised by the 
required amount. 

We shall only calculate luminosities for stars in the main sequence 
(i.e. hydrogen burning) stage of evolution. Other stages (apart from 
the final white dwarf stage) are far too complicated for the simple 
methods if calculation used here - e.g. the assumption of a simple 
density and pressure profile can hardly be expected to remain valid 
for stars with very mixed internal constitutions with heavy elements 
in the core and light ones nearer the surface. Indeed the most advanced 
stages, particularly where dynamic processes are involved, are only 
just beginning to be understood at all. However due to the relative 
rapidity of the processes involved, a star will only spend a small fraction 
of its life in these stages, so that as far as the calculations of the 
timescales of evolution in the next section are concerned, nothing will 
be lost by ignoring them. The final white dwarf stage may be comparable 
with the total evolutionary timescale but the luminosity is by no means 
steady during that period - in fact it decays by many orders of magnitude 
as the star cools. The calculations are not difficult, but we shall not 
include them here. 
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We shall start by calculating the luminosities of main sequence 
stars on the assumption that the energy is transported by radiation. 
This will give sound results so long as convection takes place only in 
restricted regions of the star. In larger main sequence stars where 
the extremely sensitive C-N-0 cycle is operative, energy will be 
generated only in a highly region near the centre and as 
a result there will be a convective core which distributes the heat over 
a large volume. Also there will in general be a thin convective layer 
near the surface of a star associated with the region where the gas 
is only partially ionised. Nevertheless, provided that there is no 
convection in the intermediate layers, the rate at which energy is lost 
will essentially be governed simply by the opacity in these intermediate 
layers. However for smaller stars with less than about | Al , 
although there is no convective core, the convective surface layer extends 
very deeply into the interior, and eventually, as was pointed out by 
Limber (1958)^^'^^ reaches the center for stars near the lower size limit. 
Under these conditions the neglect of convective transfer leads to gross 
underestimation of the luminosity so we shall give the appropriately 
modified calculation after completing the radiative case. 

The radiative energy flux vector F is given by F = (l/KnN)VPr 
where k is the opacity per nucleon per unit volume (not the 
more conventional opacity per unit mass per unit volume) and Pj. is 
the radiation pressure, i.e. by (48), ^ ^T'^ . The 
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luminosity L (i.e. the total rate of energy loss) will be given by 
L E?F where F is a typical value of the flux magnitude 
in the interior, and we may estimate the order of magnitude of 
VPr as Pr/ R where is the central radiation 
pressure. Thus we obtain 

L ^ RP, . 

In a star in which radiation pressure dominates (i.e. 77 <C 1) 
we have Pr ~ P, and so by (39) and (43) we obtain 

L ^ ^— N/N^ 

while in a star for which gas pressure dominates (i.e. 77 3> 1) 
we have, by (52), | n^^^{N/Ni^f/^ , and therefore by 

(39) we obtain 

In the transition region, where 77 1 , which occurs where ~ 10 A^l ) 
the two families agree and give L ^ 10 /{nm^) . 

It is now necessary to calculate k , the opacity per nucleon. 
In the range between the ionization temperature Ti and the pair 
creation temperature Tp , the dominant contribution at sufficiently low 
densities arises simply from Thompson scattering of photons by electrons 
for which^^ the cross section is roughly the square of the classical electron 
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radius. Therefore setting ~ as usual, we obtain 




The appearance of the mass in the dominator shows why scattering 
by heavier particles is negligible. However the nuclei do give rise 
to an important effect of a different kind when the matter is denser, 
since their electric fields will then have a significant effect on the 
electron orbits. The analysis of this situation is rather complicated, 
but approximate formulae due to Kramers are in standard use for the 
main contributions to the opacity, i.e. for those arising from bound 
free and free - free electronic orbit transitions. The Kramer free - free 
formula gives 

-1/2 

while the Kramers bound - free formula differs from this only by an 

additional factor, about 10 Z , where Z is the proportion 

by weight of elements heavier than helium, and the factor 10^ is of 

purely arithmetic origin. We have already restricted our attention 

to stars in the early stages of evolution where hydrogen is predominant, 

and in such stars Z is never more than of the order of 10~^ . 

Therefore we can use (65) to obtain the correct order of magnitude of 

the opacity even when the bound - free contribution is not less important 

than the free - free contribution. 

Therefore the Thompson scattering formula (64) should be 
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replaced by (66) whenever the latter gives a higher opacity. 

using (49) we can rewrite (66) in the convenient form (where the coefficient has 

been doubled to allow for a comparable bound - free 

contribution) 



The formula (66) differs from the formula (64) by the factor 

(T/me)~"^/^e^?7 . it will be noticed that the factor (T/mc)~^/^e^ 

is never greater than unity in the temperature range under consideration 

(it is in fact of order unity at the minimum temperature Tj e^mg 

and is less than this at higher temperatures). Therefore the formula 

(66) will always give a smaller contribution than (64) when rj 

is less than unity, or in other words electron scattering gives the main 

contribution to the opacity whenever radiation pressure is dominant. Thus when 

7] <C 1 , i.e. when is greater than about IOA^l , we shall have by (62) and (64) 



(It can be shown that this formula remains valid even when there is a 
convective core.^^) When r/ <C 1 , i.e. when is less than 
about 10 Ll , there will be two possibilities according to whether 
Thompson scattering or Kramers scattering is most important. In the 
former case we shall have by (63) and (64) 




(66) 




(67) 




(68) 
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and in the latter case we shall have by (63), (66) and (53) 

L^IQ-U—] — e-^m^'iN/N^)\ (69) 

The correct formula is whichever one gives the lower luminosity 
corresponding to the higher opacity). We see that when electron 
scattering is dominant, the luminosity depends only on the nucleon 
number, but that when the Kramer's scattering is dominant the luminosity 
depends also (but not very sensitively) on the temperature. The 
temperature we are most interested in is the hydrogen burning temperature, 
Th , at which the star spends most of its life. By (49), (69) gives in 
this case 

L ^ 10-^ (—] ' f m^-i {N/N^f . (70) 

As a consequence of (6) we see that the transition from (68) 

to (70) occurs when lOe Al ■, or, using the numerical 

value for e , N ^ Ni^ . Thus for a hydrogen burning star (67) 

gives the opacity for A^ larger than about 10 Ll , (68) gives 

the opacity for A^ in the range approximately from 10 A^l to A^l , 

and (70) gives the opacity for A^ below about Al • For the 

Sun, A^ lies just below Al as we have already remarked, and more 

accurate calculations verify^^ that the Sun does lie in the regions 

where Kramers scattering is dominant (in fact not quite so near the 

edge of this regions as these rough estimates indicate). 
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So far we have not considered the surface layers of the star, 
but it is now necessary to take them into account. The effective 
temperature, Tg , at the surface is determined by the necessity 
that the energy flux arriving at the surface should be radiated away 
into space. The energy flux F from the surface of a black 
body at temperature is given^^ by F = {'k'^/60)T^ : 
thus the effective temperature is related to the luminosity by 

Using (39) we obtain 

Substituting from (52) or (54) and from (67), (68) or 
(69) gives 

when is greater than about 10 A'^l , and 

T/ ^ lO-^e-^ ml mN (N/N^) 

or 

1 /2 

T,' ^ 10-\-' ml mN (N/N^)' (—) T' 

when is less than about 10 A^l , according to whether Thompson 
or Kramers scattering is dominant. We can use these formulae 
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to check the physical plausibihty of the picture of purely radiative 
energy transport with which we have been working. The idea that 
the star has a more or less well defined boundary, at a temperature 
from which radiation takes place into space, is perfectly reasonable 
so long as the matter of the star is effectively opaque right up to 
the boundary. However gaseous matter at low temperature has a very 
low opacity at most frequencies (except for molecular absorption lines), 
and the opacity only becomes large at temperatures where ionisation 
begins to take place. Quite suddenly, at a temperature Tq of order 
i.e. by (28) , 

To ^ lO-^Tc , 

the opacity becomes extremely large due to molecular dissociation 
and bound-bound and bound-free electronic orbit transitions. 
(At higher temperatures the opacity slowly decreases again as 
free-free transitions become relatively more important, and above 
Tj ^ |e^me it declines steadily according to the Kramers law (65), 
finally leveling off when only Thompson scattering remains important.) 

Bearing this in mind we see that it is virtually impossible for a 
star with a hot interior to have an effective radiating temperature Te 
much below To , since as the surface layers must be nearly transparent 
at lower temperatures, radiation from the hotter regions inside will 
inevitably escape, heating the outer layers in the process. Thus, if any 
of the formulae (73), (74) or (75) predicts a temperature lower 
than To , we shall have a contradictory situation. 

What will happen if this predicament arises? It appears 
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that under these conditions the star will automatically adjust itself 
so that convection takes place to a sufficient extent to raise the 
rate of energy transport to the level required to prevent the surface 
temperature dropping below Tq . The way this occurs may be 
briefly described as follows (see Stein^'^ and Speigel^^). According to 
standard convection theory, convection can take place when the actual 
temperature gradient is less than the adiabatic temperature gradient 
of the gas, i.e. when (1 — '-f~^){T/P)'VP > VT where 7 is the 
adiabatic index. The rate of convective energy transport increases 
rapidly with the excess of the adiabatic over the actual temperature 
gradient. Although the rate is very difficult to calculate, principally 
due to uncertainties in the convective mixing length, it can easily 
be shown that under the conditions obtaining in stars an extremely 
small excess temperature compared with the actual temperature 
gradient can produce a very large convective energy flux. Now in a 
star with surface temperature in the neighbourhood of To there will 
nearly always be efficient convection just below the surface because 
due to the effect of the molecular dissociation and ionisation, 7 will 
be very close to unity. If energy is not coming up fast enough from 
further inside, the outer layers will cool to temperatures not much 
greater than Tq , and the resulting reduction in temperature 
gradient will enable the convective layer to eat back further into the 
star. The final result will be that the convection will extend just 
sufficiently far into the star (if necessary going right to the centre) , 
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and will take place at just a sufficient rate, to deliver enough heat 
to maintain the surface temperature about To . The whole process 
is rather subtle, and it is only recently, largely due to the work of 
Hayashi (1962)^^, that it has come to be properly understood. 

Thus the formulae (73), (74) and (75) will only be 
correct when they predict a value of greater than Tq . 
Otherwise they must be replaced by the simple equation 

e ~ o • 

Checking this condition for the main sequence with T ^ Th ~ lO^^e 
for the central temperature, we obtain 



{Tjny^lOe-' 
when is greater than about 10 A^l , 



'2 \ 2 

e 



'2 \ 2 



when is less than about 10 A^'l , but greater than A'l , and 



(TJT,)' « lO-'e- (^) ( ^ ) {NINO' 



when A^ is below A'^l • We now notice a remarkable consequence 
of the coincidence (12) which, in conjunction with the coincidence 
(6), has the result that all the potentially large factors cancel 
out, leading to effective surface temperatures remarkably close to To 
Thus by (77) the upper limit of the surface temperature for large 
main sequence stars, i.e. blue giants, is To ^ (40 e~^)^/'^To ~ IOTq . 
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In the intermediate range of masses we have similarly 

Te ^ (4e-2) (at/aTl) 1/47; ~ 6(A^/iVL)^/^To , and finally 

for masses below JVl we have Te f« | e~^{N / Ni))^I^To ^ h{N / N^^^'^T^ 

Therefore at a critical value, N^^Ni^ the surface 

temperature drops to Tq . Stars of this order and smaller will 

transfer energy primarily by convection. By (72) and (52) we 

obtain their luminosity as 

L ^ 10-9 e'' (I^yf—V' {N/N^f 

in which form the formula applies also to larger stars with up 
to about 10 A'^L in the Hayashi convective phase as they approach 
the main sequence. On the main sequence we have T Th , so 
that for red dwarfs we can write more specifically 



lO-^e^^r^Vf^V {N/N, 



VmN/ \ rric 

The fact that the luminosities (77), (78) and (79) are 
comparable with this is only true as a result of the apparent 
fiuke that the nucleon mass is comparable with the 9th power of the 
fine structure constant. Had it been the 11th power, say, all main 
sequence stars would be convective red dwarfs. 

The results of this section are plotted as mass luminosity and 
temperature luminosity diagrams in fig 2. As in the previous diagram 
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no attempt has been made to smooth out the comers where the different 
regimes overlap, in order that the basic structure should stand out 
clearly. The two sets of results could be combined to give a 
temperature luminosity diagram - an order of magnitude explanation of 
the famous empirical Hertzsprung-Russel diagram. However although this 
diagram would be correct in order of magnitude, it would not be very 
satisfactory, because due to the remarkably small temperature range 
involved, as a result of the coincidence (12), many of the 
interesting features of the Hertzsprung-Russel diagram involve quite 
small temperature changes, which one cannot hope to account for by 
the rough methods used here. 

Having thus seen roughly how the stationary state of a star 
depends on the fundamental parameters, and knowing from earlier 
sections roughly the amount of nuclear energy available, one is 
in a position to estimate the orders of magnitude of lifetimes of 
stellar evolution. However the lifetimes, particularly of smaller 
stars, are so long that cosmological assumptions concerning the 
constancy or otherwise of the microphysical parameters are involved, and 
therefore we shall postpone this calculation to Part II of this survey, 
where more controversial matters will be discussed. 
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Postscript, October 2007. 

As the preceeding notes were too long for the journals (such as Nature) 
that seemed suitable from the point of view of subject metter, I prepared an 
abbreviated version that was ultimately published in 1973 pQ. 

It had been my intention that this "Part I", dealing with the numbers 
characterising local cosmogony, would be followed by a separate "Part 11" , dealing 
with the numbers characterising lage scale cosmology. Instead, however, the two 
parts were finally merged in a relatively short (13 page) set of lecture notes 
entitled "Large numbers in astrophysics and cosmology" that was presented at 
a Princeton meeting (organised by John Wheeler for the Clfford Centennial) on 
21 February, 1970, developing the (for some purposes indispensible [2]) notions 
that - in the version I subsequently published OH] - were designated as (strong 
and weak versions of) the anthropic principle. 

The original 1967 notes - namely the "Part I" reproduced here - were 
circulated in crude stencil printed form (c/ figures) just before the observational 
confirmation (following the discovery of pulsars) that neutron stars actually do 
exist. The epoch making introduction, at about the same time, of the term "black 
hole" was helpful for the presentation of the simpler published version [1] . The 
later did however (for the sake of wider readability) omit many of the details 
that were thought noteworthy by subsequent authors. It is for this reason, as 
well as for its purely historical interest, that the original version reproduced here 
has been directly cited in a variety of relatively ancient [3 El [71 El E] , and also 
more recent [IDl [HI [121 lEl [HI [ISl [IE] publications. The purpose of this belated 
transcription is therefore to make the omitted details more generally accessible 
via electronic archiving. It is pertinent to add some comments about how much 
of the contents of this restoration remains effectively applicable today. 

Even though it is forty years since they were first written, the main body of 
the numbered equations in this paper remains operational as a collection of rough 



order of magnitude relationships (which is all they were meant to be) holding in 
the limited regimes of application for which they were intended: only a few of 
the explanatory comments (particularly those concerning weak interactions and 
specially P processes) are now obsolete (in view of subsequent discoveries such 
as that of neutrino masses and of the W and Z bosons). 

Within the original framework, as essentially concerned only with gravita- 
tional, electromagnetic and strong interactions, there are just one or two points 
that needed rectification at the time. One item was the omission of references to 
noteworthy contributions by workers such as Ed Salpeter, and particularly Bob 
Dicke, with whom I was not yet acquainted. Another item was the over simplified 
description of the neutron drip process on page 33, where it was suggested that 
neutrons start to leak out of the nuclei as soon as the density threshold 
is exceeded, whereas this process does of course require a much higher threshold: 
its actual value is given roughly by n-^ ~ {m^/2m^)^ . 

The main caveat concerning the subsequent discussion (foreshadowing the 
strong anthropic principle) is that in the derivation of (58) it is the second ver- 
sion of (19) that should have been used, namely Th ~ 10~^e^mN . The published 
treatment [H [3] dealt with this issue correctly, but in the original treatment re- 
produced here I had naively supposed it to be the height of the classical Coulomb 
barrier, at a radius given roughly by the Compton wavelength, , of the pion, 
that controles the rate of thermonuclear reactions. It was soon drawn to my 
attention in the ensuing discussions (with experts such as John Wheeler, and 
later Paul Davies ^S\) that the rate is actually controlled by quantum pene- 
tration of the Coulomb barrier at a larger radius determined (as an application 
of Heisenberg's uncertainty principle) by the protonic analogue, e~'^'m~^ , of the 
ordinary (electronic) Bohr radius. For a properly logical treatment of this ef- 
fect, the energy factor e^m^r needs to be replaced by the numerically comparable 



energy factor 20e^mN in the equation (58), where it gives 

iVH~10-^(e'mN/me)=^/^iVL, 



(50a) 



and similarly in subsequent equations such as (70), where it gives 



and (78) where it gi 




(70a) 



(Tc/To) 



(78a) 



At the level of the weak anthropic principle this replacement makes no difference, 
because ^ 20 e^m^ in our part of the universe, but it does matter for the 
strong anthropic principle, which envisages that masses and coupling constants 
might deviate from their familiar values in other parts of what has now come 
to be known (though not unanimously loved) as a "multiverse" [19]. This rec- 
tification means that the simple expression A'^h ~ 10~^iVL for the lower limit of 
the main sequence, and the occurrence near the Landau number A^l itself of the 
transition from Thompson to Kramers opacity, are due not, as was suggested, to 
the coincidence (6) involving nuclear physical quantities, but rather to another 
coincidental relationship involving only atomic physical quantities, namely 



The (anthropic?) condition [3] that gravity be marginally weak enough for to 
drop to the order of To when A^ becomes small compared with A^l (so that heat 
transport in smaller stars will be primarily convective) can be seen from (78a) to 
be attributable, not directly to (12), but to a relation that is equivalent wherever 
- as in our own part of the universe - one has rric/m-^ ^ , namely 




(6a) 




(12a) 
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